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ABSTRACT 

Since the opacity of a magnetized plasma depends on polarization of radiation, 
the radiation emergent from atmospheres of neutron stars with strong magnetic fields 
is expected to be strongly polarized. The degree of linear polarization, typically 
~ 10 — 30%, depends on photon energy, effective temperature and magnetic field. The 
spectrum of polarization is more sensitive to the magnetic field than the spectrum 
of intensity. Both the degree of polarization and the position angle vary with the 
neutron star rotation period so that the shape of polarization pulse profiles depends on 
the orientation of the rotational and magnetic axes. Moreover, as the polarization is 
substantially modified by the general relativistic effects, observations of polarization of 
X-ray radiation from isolated neutron stars provide a new method for evaluating the 
mass-to-radius ratio of these objects, which is particularly important for elucidating 
the properties of the superdense matter in the neutron star interiors. 



Subject headings: polarization — pulsars: general — stars: magnetic fields — stars: 
neutron — X-rays: stars 



1. Introduction 

Optical and radio polarimetry has proven to be a powerful tool to elucidate properties of 
various astrophysical objects. For instance, virtually all our knowledge about the orientations of 
the magnetic and rotation axes of radio pulsars comes from analyzing the swing of polarization 
position angle within the pulse (see, e.g., Manchester &; Taylor 1977; Lyne & Manchester 1988). On 
the other hand. X-ray polarimetry has remained an underdeveloped field of astrophysics. Although 
various X-ray polarimeters have been designed (e.g., Kaaret et al. 1990; Weisskopf et al. 1994; 
Eisner et al. 1997; Marshall et al. 1998), and importance of X-ray polarimetry convincingly 



- 2 - 



demonstrated (Meszaros et al. 1988), most recent measurements of X-ray polarization has been 
made as long ago as in 1977, with the OSO 8 mission (Weisskopf et al. 1978). Nevertheless, it is 
expected that X-ray polarimeters will be launched in near future (see, e.g., Tomsick et al. 1997). 
To develop efficient observational programs for forthcoming X-ray missions whose objectives will 
include X-ray polarimetry, the problem of polarization of various X-ray sources should be carefully 
analyzed, with the main emphasis on new astrophysical information to be inferred from such 
observations. 

In the present paper we consider polarization of thermal X-ray radiation from isolated 
neutron stars (NSs) with strong magnetic fields. Recent observations with the ROSAT and ASCA 
missions have shown that several such objects arc sufficiently bright for polarimetric observations 
— e.g., the radio pulsars PSR 0833-45 and PSR 0656+14 (see Ogelman 1995, and Becker & 
Triimper 1997, for reviews), and the radio-quiet NSs RX J0822-4300 (Zavlin, Triimper, &; Pavlov 
1999) and RX J1856.5-3754 (Walter, Wolk, k Neuhauser 1996). Their soft X-ray radiation was 
interpreted as emitted from NS surface layers (atmospheres) with effective temperatures T^s in 
the range of (0.3 — 3) x 10^ K. Since photons of different energies escape from different depths 
of the NS atmosphere with temperature growing inward, the spectrum of the thermal radiation 
may substantially deviate from the blackbody spectrum (Pavlov &: Shibanov 1978). Moreover, if 
there is a strong magnetic field in the NS atmosphere, such that the electron cyclotron energy 
Ebc = heB/nieC = 11.6(i?/10^^ G) keV is comparable to or exceeds the photon energy E, then 
the radiation propagates as two normal modes (NMs) with different (approximately orthogonal) 
polarizations and opacities (Gnedin Sz Pavlov 1974). For typical magnetic fields, B ~ IQ-'-^-lO^^ G, 
the NMs at soft X-ray energies, E <^ Ebci are linearly polarized in a broad range of wavcvector 
directions, and the opacity Kg of the so-called extraordinary mode (polarized perpendicular to B) 
is much smaller than that of the ordinary mode. Kg ~ {E / Esef'i'io- As a result, the extraordinary 
mode escapes from deeper and hotter layers, so that the emergent radiation acquires strong linear 
polarization perpendicular to the local magnetic field (Pavlov & Shibanov 1978). Polarization 
of the observed radiation depends on the distribution of magnetic field and temperature over 
the visible NS surface. If these distributions are axisymmetric, the polarization is a function of 
the angle O between the symmetry (magnetic) axis and the line of sight. If the direction of the 
magnetic axis varies due to NS rotation, the polarization patterns show pulsations with the period 
of rotation, so that measuring the polarization pulse profile allows one to constrain the orientations 
of the axes. Due to the gravitational bending of the photon trajectories, the visible fraction of 
the NS surface grows with increasing the NS mass-to-radius ratio, M/R, which reduces the net 
polarization because the observer sees additional regions with differently directed magnetic fields. 
On the other hand, the gravitational field affects the magnetic field geometry making the field 
more tangential (Ginzburg &; Ozernoy 1965), which increases the observed polarization. These 
GR effects enable one, in principle, to constrain M/R by measuring the X-ray polarization. We 
demonstrate that the expected X-ray polarization of the thermal NS radiation is high enough to 
be measured with soft-X-ray polarimeters in a modest exposure time, and these measurements 
can provide important new information on both the geometry of the magnetic field and the 
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mass-toradius ratio. 



2. Description of calculations 

The intensity / and the Stokes parameters Q and C/ at a given point of the NS surface can be 
expressed as (Gnedin &; Pavlov 1974) 

I = lo + Ie, (1) 
Q = {lo -le) PL cos 2Xo , (2) 
U = {Io-Ie)PLsm2xo , (3) 

where /„ and /g are the intensities of the ordinary and extraordinary modes, pL = (1 — + 
is the degree of linear polarization of the NMs {V is the ellipticity, i.e., the ratio of the minor axis 
to the major axis of the polarization ellipse), and Xo is the angle between the major axis of the 
polarization ellipse of the ordinary mode and the x axis of a reference frame in which the Stokes 
parameters are defined. 

We calculate the local NM intensities with the aid of NS atmosphere models (e.g., Pavlov et 
al. 1994, 1995). In the present work we assume that the surface temperature is high enough for 
the atmospheric matter to be completely ionized. If the NS surface is covered with a hydrogen 
atmosphere, this assumption is justified at Tgff > 10^ K, for typical magnetic fields of NS (Shibanov 
et al. 1993). The local intensities Ig and Ig depend on the photon energy, magnetic field, and 
direction of emission. 

In the dipole approximation, valid at photon and particle energies much lower than nieC^, the 
degree of linear polarization of NMs can be expressed as 

PL = I (4) 
4 cos2 9' + g2 sin"* 9' 



where 9' is the angle between the magnetic field B and the unit wavevector k at the NS surface. 

The (angle-independent) parameter q is determined by the components of the Hermitian part 
of the polarizability tensor in the coordinate frame with the polar axis along the magnetic field 
(Gnedin & Pavlov 1974). The parameter q depends on photon energy and magnetic field (e.g., 
Pavlov, Shibanov, k. Yakovlev 1980; Bulik & Pavlov 1996). For instance, if the hydrogen plasma 
is completely ionized, and the electron-positron vacuum polarization by the magnetic field can be 
neglected, this parameter equals 

T7l/2/'T7l2 I 7T'2 TTI^ I7'„ \ T7l2 

_ {^Be + ^Bi - ^Be^Bi) - ^Be^Bi 



E'^EBe - Ebi) 



(5) 



where E' is the photon energy as measured at the NS surface, Ebi = {me/mp)EBe = 
6.32 (5/10^2 G) eV is the ion (proton) cyclotron energy. If the photon energy is much greater 
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than the ion cyclotron energy, the q parameter is particularly simple: q = Ese/E'. This means 
that the NMs are linearly polarized, pL — 1, in a broad range of directions, sin^ 0' ^ 2E'/EBe, 
at photon energies much lower than the electron cyclotron energy. It should be mentioned that 
equations (0)~(0) imply that the NM polarizations are orthogonal to each other (in particular, 
Xe = Xo i tt/2). This condition is fulfilled in a broad domain of photon energies and directions, 
except for a few special values of 6',E' (e.g., Pavlov & Shibanov 1979; Bulik & Pavlov 1996). 
Within the same approximations, the angle Xo coincides with the azimuthal angle of the magnetic 
field in a reference frame whose polar axis is parallel to k . 

To find the observed flux Fj and the observed Stokes parameters Fq and Fjj, one should sum 
contributions from all the elements of the visible NS surface. We assume the magnetic field and 
the temperature distribution are axially symmetric and define Fq and Fu in the reference frame 
such that the axis of symmetry m lies in the xz plane, the z axis is directed along the line of sight. 
In such a frame Fjj = 0, and Fj, Fq are functions of the angle between m and z. The ratio 
Pl = —Fq/Fj gives the observed degree of linear polarization, \Pl\, and the observed position 
angle: the polarization direction is perpendicular or parallel to the projection of m onto the sky 
plane for Fq > or Fq < 0, respectively. 

Since the NS radius R is comparable with the gravitational (Schwarzschild) radius 
Rg = 2GM/c^, the photon energy and the wavevector and polarization directions change in the 
course of propagation in the strong gravitational field. We will assume that the NS gravitational 
field is described by the exterior Schwarzschild solution. Since strong magnetic fields {B ^ lO^'^ G) 
are needed to obtain measurable polarization in the X-ray range, and all observed NSs with 
strong magnetic fields are not very fast rotators {P > 10 ms), the effects of rotation on the 
metric and on the observed radiation are very small. For the Schwarzshild metric, the observed 
energy is redshifted as E = QrE' , where gr = — Rg/R)^^"^ is the redshift factor. The observed 
wavevector is inclined to the emitted wavevector by the angle K — "d, k -z = cos{K — '&), where 
K is the colatitude of the emitting point in the reference frame with the origin at the NS center 
and the z axis directed towards the observer, i) is the angle between the normal to the surface 
and the wavevector direction k at the emitting point. The angle would coincide with K in flat 
space-time. In the Schwarzschild geometry, K always exceeds ■!?, i.e., some part of the NS back 
hemisphere is visible. For instance, 

j-Rg/R (Ja- 

for < vr {gr > 0.65), where a = Rsini)/ (Rggr) is the impact parameter in units of Rg (see, e.g., 
Zavlin, Shibanov, & Pavlov 1995). In particular, 



?? 1 2 

K — v ntan — I u 

2 16 



15T?-sini? 'd 

— ^ ^ + 7tan - 

sin2 7? 2 



(7) 



for u = Rg/R « 1. 
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The bending of the photon trajectories is associated with changing the direction of hnear 
polarization. The polarization direction rotates in such a way that it keeps fixed orientation with 
respect to the normal to the trajectory plane (e.g., Pineault 1977), remaining perpendicular to 
the wavevector. Without this rotation, the angles Xo and Xe would be conserved: = 4> at the 
observation point, where (f) = tan'^ {By /Bx) is the azimuthal angle of the magnetic field at the 
emitting point in the x, y, z frame. To find x"^^ with allowance for the GR effects, it is convenient 
to introduce the frame x',y', z' such that the z' axis is parallel to k , and the photon trajectory 
is in the x'z' plane (sec Fig. 1). The unit vectors along the axes of the two frames are connected 
with each other as follows 

x' = cos if cos{K — ■&) X + smip cos{K — •d)y — sin(i^ — -&) z , (8) 
y' = —siiLip X + cos ip y , (9) 
z' = cos (p sm{K — 1?) X + smip sm{K — y + cos{K — ■&)&, (10) 

where (p is the azimuthal angle of the emitting point in the x, y, z frame. Using the conditions that 
the angle between y' and the polarization direction is conserved, and the polarization direction is 

perpendicular to z at the observation point, we obtain x^^ = 9^ + where (/)' is the azimuthal 
angle of the magnetic field in the x',y',z' frame: (p' = tan~^ (By/ /Bx'). The angle (p' depends on 
K — -& and ip, and it tends to (p — <p when Rg/R — > 0. 

With allowance for the above-described gravitational effects, the observed flux Fi and the 
Stokes parameter Fq are given by the following integrals over the visible NS surface (see Zavlin et 
al. 1995): 

Fi{E,Q) = —g^J^^dfiJ^ dip{Io + Ie) = Fo{E,@) + Fe{E,e) , (11) 
FQ{E,e) = j^J dip{Io-Ie)PLCOs2{p + ct>') , (12) 

where d is the distance, /x = cos'!?, and the integrands are taken at the photon energy E' = E/g^. 

To calculate the integrands, we should know the magnetic field at the NS surface as a function 
of ji and ip. We consider a dipole magnetic field in the Schwarzschild metric. According to 
Ginzburg k, Ozernoy (1965), the field equals 

R R (2 + /)(r--m)f-/m 

^ = 2 ' ' 

where Bp is the field strength at the magnetic pole, r is the unit radius vector of a surface point, 
and m is the unit vector of the magnetic moment. The parameter 

^ , n^-2n-2(l-tx)ln(l-tx) 

accounts for the GR effect. For u < 1 (i? » Rg), we have f{u) ~ 1 + u/4 + lln^/SO. The radial 
and tangential components of the magnetic field are B^ = Bp cos 7 and B^ = {Bp/2) f sin where 
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cos 7 = r-rh = sin Q sin K cos (p + cos G cos K. Since / > 1, the GR effect makes the magnetic field 
more tangential. The projections of B onto the x, y, z and x', y' , z' axes can be easily found with 
the aid of equations 

r = sin K cos if X + sin K simp y + cos K z , (15) 
rh = sin a; + cos © z , (16) 

and equations (|8|)-(|To|). The strength of the magnetic field is 

r 1 1/2 r n —1/2 

B = {Bp/2) [(4 - f) cos2 7 + /y = Bpf [4 - (4 - f) cos^ Ob] , (17) 

where 6b is the angle between B and the normal to the surface f. 

The integration over the NS surface (eqs. ||ll|, [l^) proceeds as follows. For each point of the 
/i, if grid, we calculate the colatitude K{i!}) from equation (^) and the components of the radius 
vector r in the x,y,z and x',y'z' frames (eqs. |jl5| and 01-[|l^). This gives us the projections 
and strength of the local magnetic field B (eqs. |l^ and [jl^), the angles cp' = tan"^ (Byi / Bx'), 
9' = cos"^ {Bz' /B), and 9b = cos^^(S-f), and the degree of NM polarization pi (eq. |Q]). To 
obtain the intensities of the extraordinary and ordinary modes of radiation emitted to the observer 
from a given surface point, one needs to know the local depth dependences of the temperature 
and density in the NS atmosphere, determined by the local values of B, 9b, and T^g. These 
dependences are obtained by interpolation within a set of the diffusion atmosphere models (Pavlov 
et al. 1995). Then, the Iq and I^, intensities are computed as described by Pavlov et al. (1994) and 
Shibanov & Zavlin (1995). Subsequent numerical integration over fi,ip gives us Fj, Fq, the degree 
of the observed linear polarization and the position angle. 



3. Results 

To demonstrate how the observed linear polarization depends on photon energy, magnetic field, 
and mass-to-radius ratio, we consider a NS covered with a hydrogen atmosphere with a uniform 
effective temperature and a dipole magnetic field. We present the results for Tgfj = 1 x 10^ K, 
Sp/(10i2 Q) ^ Q 3^ i 3 and 10.0. We choose a standard NS radius ii = 10 km and three NS 
masses, M/Mq = 0.66, 1.40 and 1.92, from an allowed domain in the M-R diagram (filled circles 
in Fig. 2). These masses correspond to the redshift parameters gr = 0.90, 0.77 and 0.66, and the 
surface gravitational accelerations g/{lQ^^ cm^ s~^) = 0.97, 2.43, and 3.89. It should be noted 
that the properties of the emitted radiation are almost independent of the g value, so that the 
gravitational effects on the observed radiation are determined mainly by the redshift parameter 
Qr, i.e., by the mass-to-radius ratio. 

The left panel of Figure 3 demonstrates the observed photon spectral fluxes Fj (eq. [|ll|) 
from a NS with the magnetic field i?p = 1 x 10^^ G at the magnetic pole and the magnetic axis 
perpendicular to the line of sight, = 90°. The flux is normalized to a distance of 1 kpc. The 
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spectra are presented for the three values of the redshift parameter g,.. To demonstrate the effect 
of the interstehar absorption, we plot the spectra for the effective hydrogen column densities 
i^H = (unabsorbed flux), 1 x 10^'' and 1 x 10^^ cm~^ (the latter two are shown for = 0.77 
only). The effect of redshift is clearly seen, as well as that of the interstellar absorption: spectral 
maxima shift from 0.15 keV at nn = to 0.6 keV at nn = 1 x 10^^ cm~^. The contributions 
from the extraordinary and ordinary modes (fluxes and Fq) to the unabsorbed spectrum 
Fi = Fe + Fo are shown in the right panel of Figure 3 for gr = 0.77. At energies around the 
maxima of the flux spectra the radiative opacity of the ordinary mode significantly exceeds that 
of the extraordinary mode. Hence, the extraordinary mode is emitted from deeper and hotter 
atmosphere layers, providing the main contribution to the total flux (see Pavlov et al. 1995 for 
details). At higher energies {E' >3~^/^i?Be) the relation between the two opacities, and the two 
NM fluxes, is reversed (Kaminker, Pavlov, & Shibanov 1982). This leads to changing the sign of 
Fq (i.e., to the jump of the polarization position angle by vr/2). 

Several examples of the dependences of Pl = —Fq/Fj on photon energy are presented in 
Figure 4, for Bp = 1 x 10^^ G and different values of the angle and the redshift parameter gr- 
in the soft X-ray range, where the thermal NS radiation is most easily observed, Pl is positive, 
i.e., the polarization direction is perpendicular to the projection of the NS magnetic axis onto 
the image plane. In this energy range the ordinary mode is emitted from superficial layers with 
lower temperature, whereas the extraordinary mode is formed in deeper and hotter layers with a 
larger temperature gradient. As a result, the ratio Fe/Fo grows with E at low energies until this 
effect is compensated by the decrease of the difference between the extraordinary and ordinary 
opacities (see the right panel of Fig. 3). At higher energies the extraordinary and ordinary fluxes 
approach each other, so that Fe/Fo decreases with increasing E reaching unity at E ^ O.SgrEBe 
(hereafter, Ese and Esi are the cyclotron energies for the magnetic field Bp). Since f« 1 at 
Ebi <^ E' <^ Ebs, and if + (j)' does not depend on E' , it follows from equations ( pT[ ) and (jl^) that 
Pl oc {Ff, — Fo)/{Fe + -Fo)) with a proportionality coefficient independent of E. This explains the 
energy dependence of Pl in Figure 4. In particular, starting from energies E ~ grEsi, Pl grows 
with E until it reaches a maximum (at ~ 1 keV for Bp = 1 x 10^^ G). At higher energies the 
polarization spectra steeply decrease with increasing E, reach zero at ~ O.SgrEse, where the 
contributions from the two NMs cancel each other, and become negative (polarization direction 
becomes parallel to the m projection) at higher energies, where the flux decreases exponentially 
at the effective temperature chosen. 

As expected, the degree of polarization grows with increasing Q from 0° to 90°. It equals zero 
at = 0° (the magnetic axis is parallel to the line of sight) because the contributions from the 
azimuthal angles ip and ip + tt/2 to Fq (eq. |]l2|) are polarized in orthogonal directions, being of 
the same magnitude {By' = 0, 0' = 0, cos 2ip = — cos 2((/? it it/2)). The polarization is maximal at 
B = 90°, when the magnetic lines are seen by the observer almost parallel to the magnetic axis on 
a substantial (central) part of the visible stellar disk. 

We see from Figure 4 that the degree of polarization decreases with increasing M/R (or 
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decreasing gr). The main reason is that the observer sees a larger fraction of the whole NS surface 
due to stronger bending of photon trajectories. As a result, the overall pattern of the magnetic 
lines on the visible NS disk becomes more nonuniform, which leads to additional "cancellation" 
of the mutually orthogonal polarizations emitted from parts of the surface with orthogonally 
oriented magnetic line projections. This effect is partly compensated by the other GR effect, more 
tangential dipole magnetic field in the Schwarzschild metric (the parameter / in eq. equals 
1.08, 1.14, and 1.22 for gr = 0.90, 0.77, and 0.66, respectively), but the effect of bending prevails. 

Figure 5 demonstrates the effect of magnetic field strength on the degree of polarization. In 
particular, Pl in the soft X-ray range grows with Bp for typical NS magnetic fields. When the 
magnetic field is relatively small. Bp <1 x 10^^ G, the fast growth of Pl is due to the increasing 
difference between the extraordinary and ordinary opacities (hence, increasing differences between 
the escape depths and between the NM intensities). Slower growth of P^ at intermediate fields 
(around Bp ~ 1 x 10^^ G) is mainly due to the decrease of the surface temperature with increasing 
B (Pavlov et al. 1995), which reduces the ordinary flux Fq and the ratio Fo/Fg. With further 
increase of magnetic field, the surface temperature ceases to decrease, and Pl{Bp) saturates in 
the soft X-ray range at -Bp > 3 x 10^^ G, until the field becomes so strong, Bp > 3 x 10^^ G, that 
grEBi becomes comparable with E, and the proton cyclotron spectral feature gets into the soft 
X-ray range. 

The proton cyclotron feature in the polarization spectrum is shown in Figure 5 for 
-Bp = 1 X 10^'^ G. The shape of the feature can be explained by the behavior of the parameter q 
and the NM intensities near the proton cyclotron resonance. According to equation (^), q grows 
with decreasing photon energy until E' reaches S^^'^Esi', then it sharply decreases, crosses zero 
in the very vicinity of proton cyclotron resonance, at E' ~ i?_Bi(l -|- 2me/mp), and tends to — oo 
(q ~ —EseE'^^/ E'^ at E' <^ Esi)- This means that pi reaches zero, the NM intensities Jg and 
lo equal each other, and the integrand of equation (^) equals zero at an energy in the very 
vicinity of the proton resonance corresponding to the local magnetic field. If the direction of 
the local magnetic field is such that cos2((^ + (f)') > (see eq. [Q), which roughly corresponds 
to the projection of the local magnetic field onto the sky plane within ±45° of the magnetic 
axis projection, then the integrand in the expression for —Fq is positive at energies around the 
resonance, so that the energy dependence of the integrand for the corresponding surface points 
looks like an "absorption line" in a positive continuum, with its minimum (zero) value at the local 
resonance energy. Integration over the area with positive cos 2((/? + (^') yields a positive contribution 
to —Fq, with an absorption line somewhat broadened, and its minimum above zero, because of 
nonuniformity of the magnetic field. On the contrary, the energy dependence of the integrand in 
the area where cos2((/? + (f)') < looks like an "emission line" on a negative continuum, with its 
maximum equal zero at the local resonance energy. The integral over this area gives a negative 
contribution to —Fq, its absolute values are minimal in the energy range which includes the local 
resonance energies. If the strengths of the local magnetic fields are different in the areas of positive 
and negative cos2((^ + (p'), the integration over the whole visible disk results in a complex feature 
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in the —Fq{E) and Pl{E) spectra. In Figure 5 this feature is most pronounced at i?p = 1 x 10^^ G 
and G = 90°. It consists of two components: the absorption component with a sharp minimum 
at ~ 28 eV corresponding to the field at the magnetic equator {B^, = Bpf /2), and the emission 
component at energies below ii^ ~ 48 eV, corresponding to the field at the magnetic pole. The 
feature is also clearly seen at the same Bp and Q = 45°, whereas only the emission component 
of the feature is seen at > 10 eV in the curves of Pl{E) for Bp = 3 x 10^^ G. If the magnetic 
field is superstrong, B ~ 10^^ — 10^^ G, as suggested for anomalous X-ray pulsars and soft gamma 
repeaters, the proton cyclotron feature gets into the soft or medium X-ray range, and its detection 
would enable one to measure directly the magnetic field strength. 

Generally, the polarization spectra are much more sensitive to the strength of magnetic field 
than the intensity spectra (Fig. 6). The main effect of B on Fj{E) spectra in the soft energy 
range is a shallow proton cyclotron absorption feature (see the spectrum for Bp = I x 10^^ G in 
Fig. 6). At S ~ 0.2 — 1.0 keV the intensity spectra at different magnetic fields typical for NSs are 
almost indistinguishable, with the main contribution coming from the extraordinary mode whose 
spectrum is almost independent of B at Esi <^ E' <^ Ese- 

If the angle a between the magnetic and rotation axes differs from zero, the projection of m 
onto the sky changes its orientation with the period of rotation. This means that the angle Q and 
the polarization position angle 5 with respect to a fixed (nonrotating) direction also oscillate with 
the rotation period P: 

cos = COS C COS a + sin Q sin a cos 27r<I> , (18) 

cos a sin C — sin a cos C cos 27r<I> , , 

tan (5 = ^ -^^ , (19) 

sma sm27r<l' 

where C, is the angle between Q, (NS rotation axis) and the line of sight, $ = t/P is the rotation 
phase, the angle 5 is counted from the projection of onto the sky. Equation (|l9| ) is applicable 
in the case of polarization perpendicular to the projection of m onto the sky plane {Pl > 0); for 
Pl < 0, the left-hand-side is replaced by cot 6. 

Figure 7 shows several characteristic dependences of -Pl($) and S{^) for E = 0.3 keV, 
Qr = 0.77, and a few sets of (,a. In the particular case of an orthogonal rotator, ( = a = 90°, we 
have = = 0, i.e., the degree of polarization oscillates between zero (at $ = 0, 0.5) and a 
maximum value, Pl = 25% at $ = 0.25, 0.75, showing two pulses per period, while the position 
angle remains constant (the polarization is oriented along the direction of the rotation axis). For 
the case ^ = 60°, a = 50°, the minimum polarization, Pl = 1% at ^ = 0, 1, corresponds to 
= C — a = 10°. The polarization pulse has two maxima per period, Pl = 25% at ^ = 0.36, 
0.64 (corresponding to = 90°) and a local minimum, Pl = 22% at <I> = 0.5 (corresponding 
to = C + o = 110°). The position angle oscillates around tt/2 (or —tt/2). For ^ = a = 45°, 
Pl has one broad maximum per period, of the same height as for the orthogonal rotator, when 
= 90° at ^ = 0.5. The position angle swings from to vr during the period, crossing 7r/2 at 
the phases of maximum polarization. Finally, for = 40°, a = 10°, the phase dependence of the 
degree of polarization is almost sinusoidal, with one maximum per period; it oscillates between 
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7% and 16%. The position angle oscillates around 7r/2 (or — 7r/2) with a small amplitude because 

of the small value of a. Thus, we see that a variety of the phase dependences of Pi and 5 can be 
obtained for various combinations of i^, a, which is potentially useful for evaluating these angles 
from polarimetric observations. 

4. Discussion 

As we see from the examples above, the degree of linear polarization of thermal NS radiation 
is quite high, up to 20%-50% in pulse peaks, and about twice lower for the phase-averaged 
polarization, for typical NS parameters. Optimal energies for observing the polarization are in 
the soft X-ray range, ^ 0.1-1 keV, for typical surface temperatures of young and middle-aged 
NSs. The polarization can be observed at these energies with the use of multilayer coated 
mirrors which provide high reflectivity at large grazing angles — see Marshall et al. (1998) for a 
concept for a satellite-borne polarimeter (PLEXAS) which would be able to measure the linear 
polarization from brightest thermally emitting pulsars to an accuracy of l%-3% during modest 
exposure times ~ 30-100 ks. The polarization can be measured in several narrow energy bands. 
Complementing the spectral flux with the polarization spectral data and comparing the both with 
the NS atmosphere models would considerably narrow allowed ranges for NS parameters such as 
the magnetic fleld and effective temperature. Particularly strong constraints could be obtained if 
the X-ray polarimetric measurements are supplemented by measuring optical-UV polarization of 
the same sources. Several middle-age pulsars have been observed successfully in the optical-UV 
range with the Hubble Space Telescope (e.g., Pavlov, Welty, Sz Cordova 1997, and references 
therein). Although the expected polarization in this range is lower than in soft X-rays, it still can 
be as high as 5%-15% (see Figures 4 and 5), so that UV-optical polarimetric observations of these 
sources seem quite feasible. 

Using detectors with high timing resolution (e.g., microchannel-based photon counters) for 
polarimetric observations of pulsars would allow one to measure phase dependences of the degree 
of polarization and the position angle. These observations would be most useful to determine the 
inclinations of the rotation and magnetic axis {( and a). Although radio band polarization data 
have been widely used to constrain these angles, the results are often ambiguous because the same 
behavior of the position angle can be fitted with different combinations of C and a. An additional 
difficulty in interpreting the radio polarization data is caused by shortness of the duty cycle of 
radio pulsars the radio flux is too low during a substantial fraction of the period to measure the 
polarization. Since the X-ray flux remains bright during the whole period, measuring the position 
angle and the degree of polarization in X-rays would enable one to infer the orientation of the 
pulsar axes with much greater certainty. 

Of particular interest is the result that the observed polarization is sensitive to the NS 
mass-to-radius ratio, the most crucial parameter to constrain still poorly known equation of state 
of the superdense matter in the NS interiors. Since the radio emission of pulsars is generated well 
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above the NS surface, this ratio cannot be constrained from polarimetric observations of pulsars 
in the radio band. Although some constraints can be obtained from the pulse profiles of the X-ray 
flux (Pavlov &; Zavlin 1997), the polarization pulse profiles are more sensitive to the gravitational 
effects. 

Primary targets for studying the X-ray polarization of thermal NS radiation are the X-ray 
brightest, thermally radiating pulsars such as PSR 0833-45 (Vela) and PSR 0656-1-14. A prototype 
of another class of promising targets, radio-quiet NSs in supernova remnants, is RX J0822-4300 
in Puppis A. X-ray polarimetric measurements would be crucial to establish the strength and 
geometry of the magnetic field of this putative isolated X-ray pulsar (Pavlov, Zavlin, k. Triimper 
1999). Also it would be very interesting to study X-ray polarization of X-ray bright, radio-quiet 
isolated NSs which are not associated with supernova remnants. Particularly interesting is 
the object RX J1856. 5-3754 which has a thermal-like soft X-ray spectrum but does not show 
pulsations (Walter et al. 1996). The lack of pulsations may be explained either by smallness of the 
magnetic inclination a, if the magnetic field is typical for NSs (~ 10^^-10^^ G), or by a very low 
surface magnetic field. Polarimetric observations would enable one to distinguish between the two 
hypotheses — the polarization is expected to be high (and unpulsed) in the former case (unless C, 
is also small), and it would be very low if the magnetic field is lower than 10"'^'' G. Distinguishing 
between the two options is needed to choose either low-field or high-field NS atmosphere models 
— applications of these types of models to interpretation of the multiwavelength observations of 
this source yield quite different NS parameters (Pavlov et al. 1996). 

It follows from our results that thermal radiation from millisecond pulsars, whose typical 
magnetic fields are 10^ — 10^ G {Ebs ~ 1 — 10 eV), is not polarized in the X-ray range. However, 
their polarization is expected to be quite strong in the optical-UV range and can be measured in 
sufficiently deep observations. The best candidate for such observations is the nearest millisecond 
pulsar J0437-4715 whose magnetic field, i? ~ 8 x 10^ G, was estimated from radio observations. 
The thermal radiation from the NS surface, which is expected to be heated up to ~ 10^ K, should 
prevail over the radiation from the very cool white dwarf companion at A < 2000 A. This means 
that polarization from this pulsar could be observed in the far-UV range with the Hubble Space 
Telescope. It should be mentioned that the rapid rotation of millisecond pulsars may affect not 
only the intensity pulse shape (Braje, Romani, &; Ranch 1999), but also the polarization, the effect 
neglected in the present paper. 

In summary, our results demonstrate that including X-ray polarimeters in future X-ray 
observatories, or launching dedicated X-ray polarimetry missions, would be of great importance 
for studying NSs. The polarimetric observations would be useful for studying not only the thermal 
component of the NS X-ray radiation, but also the nonthermal component which dominates in 
many X-ray emitting pulsars, particularly at higher energies. Further theoretical investigation 
of polarization of both thermal and nonthermal X-ray emission from NSs are also warranted 
to provide firm interpretation of future observational data. In particular, it would be useful to 
consider the polarization with allowance for possible nonuniformity of the temperature distribution 
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over the NS surface and to study the effects of chemical composition of the NS atmosphere on 
polarization. 

We are grateful to Hermann Marshall for the useful discussions of capabilities of modern soft 
X-ray polarimeters. This work has been partially supported through NASA grant NAGS- 7017. 
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Fig. 1. — Angles and vectors used in the paper. 

Fig. 2. — Mass-radius diagram for NSs. The thick sohd curves show the M{R) dependences 
for several equations of state of the superdence matter (Shapiro &; Teukolsky 1983): soft (tt), 
intermediate (FP) and stiff (TI and MF). The thick straight line, R = 1.5Rg, corresponds to the 
most conservative lower limit on NS radius at a given mass. Thin dashed lines correspond to 
different values of the redshift parameter Qr (the numbers near the lines). The filled circles give 
the NS mass and radius used in our computations. 

Fig. 3.— Left: Photon spectral fluxes Fj{E) for d = 1 kpc, Bp = I x 10^^ Q, & = 90°, and three 
values of the redshift parameter, gr = 0.66, 0.77, and 0.90 (dash-dot, solid, and dashed curves, 
respectively). The numbers near the curves denote the interstellar hydrogen column density uh, 
in units of 10^° cm~^. Right: Contributions from the extraordinary (e) and ordinary (o) NMs to 
the observed flux for gr = 0.77, nn = 0. 

Fig. 4. — Effect of the mass-to-radius ratio and the inclination O of the magnetic axis on the 
polarization spectrum for Bp = Ix 10^^ G. Thin, medium, and thick curves are for g^ = 0.90, 0.77, 
and 0.66, respectively. 

Fig. 5. — Effect of the strength of magnetic field on the polarization spectrum for g,. = 0.77. The 
curves are plotted for = 90° (solid) and 45° (dashed), and Bp/{10^^ G) = 0.3, 1.0, 3.0, and 10.0 
(the numbers near the curves). 

Fig. 6. — Unabsorbed spectral fluxes Fj{E) for gr = 0.77 and Bp/10^^ G = 0.3 (short dashes), 1.0 
(solid), 3.0 (long dashes), and 10.0 (dash-dots). 

Fig. 7. — Dependences of the degree of polarization and the position angle on NS rotation phase ^ 
for gr = 0.77, E = 0.3 keV, and different angles between the rotaiton axis and the line of sight, and 
between the magnetic and rotation axes: C, = a = 90° (solid), ( = a = 45° (long dashes), C = 60° 
and a = 50° (dot-dashes), and C = 40° and a = 10° (short dashes). 
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